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Introduce ion 


In  this  article  we  shall  consider  random  variables  (r.v.'s)  which  take 
values  on  the  closed  interval  [0,"].  We  shall  regard  these  as  the  service 
life  of  some  mechanism  according  to  a  predetermined  criterion  or  the  time  of 
death  of  a  biological  entity  and  we  shall  call  these  r.v.s  lije  lengtns. 

Letting  X  denote  such  a  random  variable  with  distribution  function  F. 
We  shall  denote  the  survival  probability  (reliability)  by  F  where 

P[X>t]  =  F (t)  =  1  -  F (t)  for  t  >  0. 

We  shall  classify  such  a  r.v.  by  the  behavior  of  its  failure  rate  (or 
hazard  rate)  which  is  defined  when  a  density  f  =  F'  exists  by 

<1  •  f/F 
t 

or  by  letting  <?  =Jo  Wc  -wVe  known 


WV  she.,  always  assume  cun-  s1--  wh*. .  e  .  i  . mpact  space  of  non 

negative  integrable  ...c  cions  and  we  a-dre*-.!*  jurativea  <-c  a  a  problem  of 
estimating  n  from  enervations  on  t.w  1-I<-  ^ength. 

no  .  _c  make  predictio:  jut  life  or  death  as  in  reliability 

stud-  or  -  acement  policie  one  is  naturally  interested  in  classifying 
the  t-e-  of  random  variables  fro...  ...—cfr  the  observations  are  obtained.  This 
claao-.  ication  can  often  be  done  more  naturally.,  by  tne  behavior  of  t  -.r 
failure  rate  rather  than  the  density.  Consequently  there  can  be  pra*.  .ical 
interest  in  estimating  it. 

iti  actuarial  studies  the  eu. —....cion  o  uilure  rates,  there  called  t..._ 
force  of  mortality,  has  been  done  for  years  because  of  the  interest  in  cal¬ 
culating  the  probability  of  death  given  survival  to  a  fixed  time.  The  famous 


paper  of  Gompertz  [6]  in  1825  proposed  an  exponential  form  for  the  failure 
rate.  This  was  modified  in  1860  by  Makeham  through  the  addition  of  another 
parameter.  Extensive  applications  of  uie  Gcmpertz-Makeham  hazard  rate  to 
human  data  have  continued  to  be  made,  see  Cramer,  Wold  [5].  More  recently 
studies  by  Jones  [8]  and  Taylor  [10]  on  the  mechanism  of  aging  in  humans  was 
carried  our  by  comparisons  of  the  failure  rate.  Thus  the  use  of  failure  rate 
as  a  means  of  classifying  data  has  a  long  and  continuing  history.  Moreover 
recent  studies  which  classify  distributions  according  to  their  failure  rates 
have  been  shown  to  have  interesting  properties.  In  this  regard  we  mention  the 
work  of  Barlow,  Marshall  and  Proschan  who  have  shown  e.g.  that  monotone  failure 
rates  generalize  some  of  the  properties  of  Polya  frequency  functions,  see  [2] 
and  [ 5  ] . 

Here  we  consider  the  problems  of  estimating  the  failure  rate  for  more 
general  classes  than  parametric,  examples  of  which  are  Gompertz-Makeham  and 
Weibull,  and  which  have  not  been  previously  studied. 

Some  of  the  useful  classifications  of  failure  rate,  wcj.1  known  by  their 
acronyms,  are  the  IF?,  (increasing  failure  rate)  and  IFRA  (increasing  failure 
rate  average).  Following  the  notation  pre/iously  introduced  in  [10]  we  .ay, 
for  any  positive  functions  q  and  r  de^nc..  us  [0,°°], 


(1.1) 

q  c  Sf  iff 

q/r 

is 

increasing 

(1.2) 

q  e  C  iff 

q/r 

is 

conve..  .creasing  and 

lim  £(> 0  .  0 
x+Q  r(x) 


Q  e  S1  iff  q  is  IFRA. 


Ocher  cases  may  arise 


when  the  failure  rate  decreases  and  this  may  occur 
in  a  variety  of  practical  situations.  We  thus  introduce  the  notation 

a  =  set  of  all  non-negative  integrable  functions  of  such  that 

JQ  q(x)dx  -  ». 

-  {q  C  nQ:q  e  +)  -  s0 

-  {q  e  flQ:q  is  convex} 

Other  subsets  may  be  of  interest  in  certain  applications. 


The  Likelihood  Function 


Let  an  ordered  sample  <  . . .  <  of  life  leugths  with  failure 

rate  q  be  given.  The  joint  density  of  the  sample  is, 

n  -Q(x.) 

TTq(x,)e 

i-1 

Taking  logarithms  and  setting  £  -  (x.  , . . .  ,x  )  we  obtain 

i  n 


(2.1) 


L*(qU)  -  j  [In  q(x.)-Q(x  )) 
i-1 


which  is  called  the  likelihood  function.  In  the  usual  sense  we  want  to 
maximize  this  with  respect  to  q  c  A^Csome  subset  of  Qq) . 

It  can  be  seen  that  th<^  right  hand  side  of  (2.1)  can  be  made  arbitrarily 
large  by  taking  q(xn)  large  but  choosing  q  so  that  Q(*n)  remains  finite. 
This  no  maximum  likelihood  estimate  of  q  exists,  in  the  usual  sense,  over 


Alternatively  we  consider  the  estimate  q,  defined  by 


(2.2) 


q  (x)  -  °°  for  x  >_  x 


and  q  defined  as  (O.x^),  from  the  almost  likelihood  equation 


(2.3) 


Uq|0  -  2  4n  q(x  )  -  £  Q(*.), 
i-1  i-1 


by  the  inequality 


L(q|0  >_L(q|£)  for  all  q  e  fl. 


Ke  shall  call  q  so  defined  the  almost  maximum  likelihood  estimate  (AiMLE) 


Remark  1:  For  a  given  sample  5,  L(*|0  is  a  concave  function  over  the 


convex  set  and  is  almost  surely  continuous  in  the  weak  topology. 

Proof:  That  is  a  convex  set  is  obvious.  Note  that  the  second  term  in 

the  definition  of  L(*|0  is  linear.  Thus  to  show  for  te  (0,1) 

L(tqi+(l-t)q2|0  >_  tL(qxU)  +  <l-t'L(q2|0 

it  is  sufficient  to  show  that 

n-1  n-1  n-1 

£  «.n[tq1(xi)+(l-t)q2(xi>]  >  t  I  in  q^x^  +  X  ln 

But  this  follows  since  the  logarithm  is  a  concave  function.  Let 
q^  -*■  q  i.e.,  this  implies  -*•  Q  and  since  the  logarithm  is  con¬ 
tinuous,  unless  some  component  of  £  falls  in  the  exceptional  set 
(which  will  happen  with  probability  zero),  then  L(qn|0  -*■  L(q|0.|| 

Since  for  any  given  £  the  function  L(*|0  is  almost  surely  a  con¬ 
tinuous  function  over  any  one  of  the  compact  subsets  °f  the 

corresponding  AMLE,  say  q^,  exists  by  the  preceeding  remark  and  it  is 
given  by 

(2.4)  L(q.  JO  ■  sup  L(q|0. 

qenk 

Let  n  *  (y^,...,yn  ]_)  b®  an  ordered  (n-l)-tiple  of  non-negative 
numbers  and  set 

(2.5)  Qk(n)  -  (qefik:  q(x1)-yi  i-J., . . .  ,n-l}  , 

which  being  a  closed  subset  of  a  compact  set  is  compact.  Let 
Ek  “  ^n:  )  *  then 


_ 


sup  L(qjs)  -  sup  sup  L(q|£) 
qeQ.  neE.  |qen.  (n) 

(2.6)  *  *  K 

n-1  n 

-  sup  £  Iny.  -  inf  £  (n-i+l)A.  (q) 

ne^  1  qenk(n)  1 


where  we  set 


Ai (q)  “Jx^  q(x)dx  for  i  -  with  Xq  *  0. 


a  3.  The  AMLE  in  the  case  of  increasing  hazard  rate 

In  this  case,  clearly,  A^q)  is  minimized  over  fi^(n)  by  setting, 
for  each  i  ■  1 . n-1 

q 00  ■  yi_1  for  x  x  <.  x  <  x1  with  yQ  -  0. 

Therefore 


A^q)  "  yi-l^Xi~*i-l^  1  "  1 . n* 

Thus  we  want  to  evaluate,  from  (2.6), 

n-1 

X  in  yt  -  (n-i)y1(xi+1-xi) 

where 

E1  ■  (<>,1 . Vl):  'li'!  i  •••  iVih 

To  do  this  we  quote  a  result  due  to  Brunk  et  al  [1]. 

If  h(a^,...,an>  is  a  bounded  strictly  concave  function  for 
-®  <  a^  <  •  i»l,...,n  which  is  to  be  maximized  over  the  convex  set 


(3.1) 


sup 

neE. 


(3.2) 


S  •  { (a. , . . . ,a  )•  a,  '  fln  ^  ...  ^  a  ] 
n  1  n  1—2—  —  n 


then  it  has  a  unique  maximum  which  we  designate  by  a  ■  (a a^).  More¬ 
over,  if  a**  -  (a^,...,ajj)  is  the  unconstrained  maximum  then  we  have  the 


following  relationship  holding:  If  a^  <  a^+^  f°r  some  k  -  l,...,n 


then 


0-  -  0  ,  -  ^  0  -  0 
a.  <  a  ,<  a,  , .  <  a,  .  and  a  >  a  ,  a.  <  a,  . 
k  —  k  k+1  —  k+1  n  —  n’  1-1 


The  usefulness  of  this  theorem  lies  in  the  following  procedure  which  is 

called  "Brunkizing"  see  [3]  and  (A),:  First  obtain  a®  the  unconstrained 

0 

maximum  of  h  over  R  .  If  a  e  S  , 

n  n 


it  is  the  solution  a  sought. 


and 


If  it  does  not  fall  in  S  then  for  some  k  ve  Duet  have  aP  >  aP  , 

n  k  k+1 

hence  by  the  theorem  above  follows  a^  ■  Thus  we  identify  the  vari¬ 

ables  a.  and  a,  , ,  to  obtain  a  new  function  of  (n-1)  variables.  We 
k  k+1 

again  maximize  this  function  over  ^ ,  checking  to  see  if  this  unconstrained 

maximum  falls  in  S  .  .  If  it  does  we  are  through.  If  not,  we  repeat  this 

n-i. 

procedure  until  at  the  final  stage  (which  must  occur  in  at  most  n  steps) 
we  have  maximized  in  unconstrained  fashion  a  function  of  a  reduced  number 
of  variables.  We  then  set  equal  to  the  final  value  of  the  pooled 
variate  to  which  that  coordinate  has  contributed. 

In  many  instances  the  unconstrained  maximum  may  be  a  stationary  point 
of  the  function  h  which  can  be  found  by  setting  the  partial  derivatives 
equal  to  zero,  i.e.,' 

Dih(a)  -  0  for  i  -  l,...,n 

and  solving  for  the  solution  -  (afj\...,a^).  The  procedure  then  is  clear. 
This  is  especially  useful  when  h  is  separable. 


Lemma:  If  a  convex  function  f  is  separable,  i.e.,  of  the  form 

n 

f(x1...,xn)  -  £ 

where  each  f^  is  convex  over  [0,«]  and  we  let  u.(i-k,i+J)  be  the 
value  if  x  which  minimizes 


i+J-1 

I 

a-i-k+1 


f  <x) 
a 


then  the  minimizing  point  of  f  over  the  region  0  <  <_  ...  <_  x^  <_  «*, 

call  it  (a^,...,a  )  is  given  by 


a  ■  max  min  u(i-k,i+J). 
k>l  J>1 


This  is  essentially  the  method  used  by  Brunk  loc.  cit.  to  find  some 
maximum  likelihood  estimates  of  parameters  constrained  in  the  same  manner 
as  (3.2). 

Moreover  applying  this  procedure  to  (3.1)  leads  to  estimates  q  which 
have  been  obtained  by  Grenander  [7]  and  studied  in  detail  by  Proschan  and 
Marshall  [9].  We  state  these  results  as 

Theorem  1:  If  Q  e  C^,  i.e.,  q  is  IFR,  then  based  on  the  sample  £  the 

AMLE  of  Q  is  Q(  •  |  £) ,  which'  is  the  greatest  convex  rainorant  of 

{t,-«n[l-F  ( t ) ] ) :  t>  0}  where  F  is  the  empirical  cumulative  dis- 
n  n 

tribution  based  on  £  and  q  is  given  explicitly  by 


5(xJo 

q(xj£) 


min  max 
j>k+l  i<k 


j  (n-ot+1)  Ax. 


t1 


o-i+i" 


j-i 


for  k  -  1, . . . ,n-l. 


It  is  quite  apparent  that  the  estimation  problems  which  arise  for  certain 
of  the  separate  subclasses,  as  defined  in  Section  1,  are  related.  We  now  state 
Theorem  2:  If  q  e  Sr,  with  r  known,  then  the  AMLE  of  Q,  based  on  £  is 

Qvt | £>  -  Q(R(t) | R(£) ) 

where  R(£)  ■  (R(x.  ),..., R(x  ))  . 

'in 


Proof:  By  definition,  and  a  slight  abuse  of  notation, 

n-1  n 

L(q|£)  -  ^  in  q (x , )  -  X 
i-1  1  1 


n-1  ,q  (x  )  v  n-1  n  rx  ,  . 

'  A  tnfc<v)  + 1 in  r<Xi>  ■  |  L 


7 


Note  that 


•L(qU)  -  L(QR~1:g(0)  “  1  In  rR_1(x. ) 

1  1 

where  juxtaposition  of  function  denotes  composition. 

Since  q  e  Sr  iff  QR  ^  e  Cq  it  follows  that 

n-l  , 

sup  L(Q|0  -  sup  L(iJ)  j  R(£)  )  -  £  £n  rR  (x.). 

{Q:qeSr>  ^eCq  1 

Hence  Q(*|£)  is  found  at  the  AMLE  <K*|R(0)  in  Cq  but  since 
^R  -  Q  we  have  0(R(t)|R(O)  ■  Q(t|0.  || 

Thus  having  an  explicit  expression  for  the  AMLE  of  q  when  it  is  known 
to  be  increasing  enables  one,  using  this  result,  to  immediately  calculate 
the  AMLE  for  q  e  for  any  known  function  r. 

As  an  illustration  we  suppose  q  e  i.e.  we  know  the  failure  rate 
divided  by  time  is  an  increasing  function.  Thus 

2 

r(x)  •  x,  R(x)  ■  for  x  >  0. 

By  the  theorem  above  we  have  our  estimate  q  defined  for  t  >  0  by 

n 

q(t|0  -  I  q[R(xk)l?<0]r(t)/xk<t<xk+1J 
k«l 

where  we  set  xn+1  *  “  and  (?(  is  the  indicator  of  the  relation  P  being 
one  if  true  and  zero  otherwise.  More  explicitly 

V1  min  max  ^  (n-a+1)  (x^-x^  .) 

lc-1  Jik+1  **  a .1+1  - 2(3-1)  K-C<Vl 

-  - 

t  >  X 

—  n 


S0 


Because  Q  e  SR  iff  QR  *  e  and  Q  e  iff  q  is  IFRA  we  see 
chat  if  ve  could  find  a  maximum  likelihood  estimate  for  q  when  it  is  IFRA 
we  could,  in  a  manner  similar  to  the  theorem  above,  find  a  corresponding 
estimate  for  Q  c  SR  for  any  known  R.  Unfortunately  we  cannot  do  so,  as 
ve  now  show,  because  the  likelihood  equation  contains,  for  each  observation 

in  q^)  -  Q(xk)  "  in^G’  (xk)+G(xk)]  -  xR  G(xk> 

where  we  define  the  increasing  fucc'tion  G(x)  -  Q(x)/x.  Clearly  we  can  make 
the  slope  G'(xk)  arbitrarily  large  while  holding  G(xR)  finite.  Thus  no 
likelihood  estimate  exists  since  this  holds  for  every  observation. 


// 


§  1  ,  The  AMLE  in  the  Case  of  Convex  Hazard  Rate 


Let  ^2  be  the  subset  of  ^  each  element  of  which  has  a  derivative 
existing  a.e.  Since  this  subset  is  dense  in  from  (2.4),  we  have  that 


(4.  )  L(q2|0  ■  sup  L(q|s). 

<1^2 

* 

Let  ft ^  t>e  the  subset  of  each  element  of  which  is  piecewise  linear, 

then 


sup  L(q|0  >_  sup  L(q|t). 

qefl^  qefl£ 


How 


(4.3) 


n-1  n 

sup  £  in  y .  -  inf  £  (n-i+l)A  (q) 

neE^  1  qefl^n)  1 


where 


E^  ”  **  ^  ’ 

Remark:  If  q  e  fl^n)  then  there  exists  a  piecewisp  linear  convex  function, 
say  p^,  defined  by 

Pq(x)  -  maxty^q' (xi)  (x-xA) :  i-l,...,n-l} 
and  such  that  for  each  i  ■  l,...,n 

(4.4)  A^(q)  >  A^p^)  with  equality  iff  q'(xi)  -  q'(xi_1). 

Proof:  The  inequality  (4.4)  is  geometrically  obvious.  If  q'(x^)  <  q'(x^) 
then  the  lines  d^x)  -  +  q'(x^)(x-x^)  and  dj_^(x)  intersect  at 

some  point  s^  and  x^  ^  <  s^  <  x^.  Now 

q'(t)  >  qMx^)  ■  Pq(xi_i)  for  xi_i  <  t  <  8a 


42  . 


and  therefore 


q(x)  -  q(xiX)  mJx  q'(t)dt  1  qp^.iHx-Xj)  -  qp(x)  -  q^*^)* 


By  similar  reasoning  we  obtain  the  Bame  inequality  on  (s^,x^).  Thus 


q(x)  >  q  (x)  for  x  e  (x.  -,x. ) 
—  p  l—i  1 


and  integration  yields  the  result.  The  claim  for  equality  is  easily 
seen.  | | 

However  from  this  remark  we  notice  that 


sup  L(q|0  £  sup  L(p|0 

q  e«2  pefi2 


and  thus  by  (4.2)  we  have 


(4.5) 


L(qJO  -  sup  L(q|U 
qefi* 


sup  sup 
neE£  pefi*(n) 


n-1 

I  An  y 


n 


2  (n-i+l)Ai(p)j 


where 


e2  ■  (n:  <n)  i  0} 


For  given  t,  *  (x.,...,x  ),  let  r)  ■  (y, » •  •  •  »y  , )  e  E-  be  such  that 
in  i  n— l 

* 

there  exists  p  e  fl^Cn),  which  is  a  piecewise,  linear,  convex  non-negative 
function  with  slope  when  passing  through  the  points  (x^y^  for 

i  ■  l,...,n-l.  Consequently  we  must  have,  as  a  necessary  and  sufficient 


condition, 

(4.6) 


n-1 


&y2  Ay, 

m,  <  i —  <  m»  <  ...  <  t - 

1  -  4x2  -  2  -  -  4Vl 


<  m  . 
—  n-1 


/3- 


Let  us  sec 


(4.7) 


xn  -  0,  y.  **  y,  -  m,x,  ,  y  -  y  ,  +  m  ,  Ax  ,  m  ■  m  . 
0  ’  J0  11*  'n  'n-1  n-1  n  n  n-1 


(4.8) 


Define  che  salient  points  of  p  as 
-Ayi  +  A(m^x^) 

s ,  ■  - - -  for  all  i  such  that  m.  >  m.  . . 

i  Am.  i  i-1 


Now 


rx.  x  rx. 

A.  (p)  mJ  P (t)dt  >*  t  p(t)  -J  t  dp(t)  , 

1  xi-l  xi_i  xi-l 

and  using  the  fact  that  p'(t)  -  m^  for  t  e  where  we  set 


s  -  x.  -  0,  s  -  x  , 
1  0  n  n 


we  have  for  i  ■  l,...,n. 


/: 


^  i  1  rf  ry  in,  n  r\ 

t  dp(t)  -  ——  (s^-x^)  +  -j-  (x^-8t)  . 


'i-1 

Substituting  into  equation  (4.5)  we  find  the  expression  in  braceB,  which  is  to 
be  maximized,  becomes 

,  2  . 
n-l  n  /  x ,  v  n-i  # , ^  a 

(4.9)  I  ^  |  8iAmi  * 

Remark:  Equation  (4.9)  defines  a  concave  function  of  the  vector 


-  (y 


*  *  •  ! 


Vi’V 


’Vi)f 


Proof:  That  the  first  summation  is  concave  and  the  second  linear  is  obvious. 
We  show  that  the  third,  it  being  subtracted,  is  convex,  by  examining 
each  term,  namely 


2. 


( y i  ^-y^+m^x^-m^  * 


ni  - 


i-1 


for  1  *  2 | i « • f n-le 


y/. 


1 

! 


To  see  this  is  convex  we  need  only  see  that  for  0  <  t  <  1  and 


l>i  >  b2  >  0  we  have 


(ta.+(l-t)a,)2  a?  a2 

tb:  +  a-t)b2  - c  bx  +  u  t}  b2 


which  can  be  easily  checked  by  elementary  methods.  | | 

We  have  just  argued  that 

L(q2|0  ■  sup  h(a) 
aeS2n-2 

it 

where  S2n_2  is  the  convex  set  of  vectors  a  satisfying  (4.6)  and 
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Lemma:  From  (4.8)  there  follows 
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Remark: 


We  now  find,  recalling  s^  ■  0,  sfl  -  xn>  that 
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Proof:  By  straightforward  calculation  we  find 
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which  simplification  shows  to  be  equa)  to  the  result  given,  and 
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from  which  the  presented  result  follows.  The  special  cases  k  ■  1  and 
k  =  n-1  can  also  be  shown  to  yield  the  result  given  by  using  the  boundary 
conditions.  |  | 

We  want  to  find  the  solution  to  the  simultaneous  equations 


V°  J-x . n-1;  V°  . n‘1- 

This  is  equivalent  with  finding  the  solution,  in  the  variables  (y^, . . . ,yn_^) 
(sj,...^  ^),  to  the  equations 


<4-10)  xj '  sj  •'/s^  (sj«'V 


J  “  1 f • • • t n— 1 


(4.11)  —  -  (n-j)  (8j^-Xj  )  +  (n-j+l)(Xj-8  )  j  - 

Is  it  possible  that  such  a  solution  exists?  The  answer  is  no  if  one  imposes 

the  conditions  s.  ■  0,  s  -  x  .  This  means  that  h  does  not  attain  a  maximum 
1  ’  n  n 

over  the  unconstrained  variables  and  so  Brunkizing  the  function  is  not  possible 
in  the  usual  sense. 


Let  n  be  fixed  and  set  X.  ■  >/  nTj,~  for  j  ■  l,...,n-l,  X  -  0. 

J  n-j+1  n 

Then  the  set  of  equations  (4.10)  is  equivalent  with  the  system 


(4.12) 


s  j  -  Xj  (1+Xj)  -  Vj8j+1  J  “  lt..., n-1. 


Introduce  the  new  variables  t^  ■  a^s^  for  1  ■  l,...,n,  where 
a^+1  «  ~aj^j  ^°r  ^  *  l,...»n-l.  By  taking  the  product  of  both  sides  for 

j  ■  we  obtain 
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where  a^  >  0  is  unspecified.  Multiplying  (4.12)  by  a^  yields 
Atj+^  ■  a^(l+Xj)x^  for  j  -  l,...,n-l.  Summing  over  j  -  l,...,k-l,  we 
obtain,  recalling  —  0 

k-1  a  (1+X  ) 

sk  -  1  „ - L-  x  k  -  2 . n-1. 

K  j-1  Sk  J 

Substituting  the  definitions  of  the  constants,  simplifying  and  rearranging 
the  summation  yields 


k-2 


(4.13) 
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On  the  other  hand  summing  over  j  ■  k,...,n-l 


the  fact  that  s  ■  x  we  find 
n  n 


k  *  2, . . . , n-1 . 
in  the  same  way  and  using 


(4,1A)  8k  "  xk  ~  2  ("1)J  ^n-kil  Axj+1  k  “  2 . n_1* 

j-k 

For  consistency  we  must  have  these  two  expressions  yield  the  same  solu¬ 
tion.  This  is  equivalent  with 


T  (-l)‘)’k  v/— fj  ■  Ax, . .  -  0  for  all  1  <  k  <  n  <  », 

jto  n_k+1  J+1  “  “ 

which  will  be  false  except  for  a  set  of  probability  zero. 

Of  course  this  means  that  no  maximum  exists  for  the  function  h  over 

the  positive  orthant.  It  must  be  that  the  maximum  will  occur  on  the  boundary 
* 

of  $2n  2  and  thus  at  some  point  we  must  consider  the  equations  with  certain 


of  the  variables  identified. 


Lee  us  suppose  chat 
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and 


>w  ■  yk +  \(\+rxk)  1  • 0 . r 


Then  writing  the  function  h  in  terms  of  these  reduced  variables  we 
find,  after  lengthy  but  straightforward  calculations  similar  to  those  done 
previously, 


(4.15) 


(4.16) 
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+  (n-k+1)  sk  -  (n-k-r)sk+r+1  . 


§5.  Some  Critical  Remarks 

It  is  now  clear  that  those  concave  likelihood  fui  ctions  which  are  to  be 
maximized  over  the  region  0  <_  x^  x^  <  «  and  which  are  not  separable 

so  that  Lemma  1  can  not  be  utilized,  and  further  do  not  possess  a  maximum 
over  the  positive  orthant,  so  that  Brunkizing  can  not  be  utilized,  pose  special 
difficulties  in  obtaining  explicit  formulae  for  the  maximizing  functions. 
Unfortunately  this  is  the  case  for  convex  failure  rates  in  general.  It  is 
clear  irom  the  preceeding  equations  that  in  each  instance  the  AMLE  estimates 
coul  1  be  found  using  a  machine  program.  This  would  begin  by  assuming  the  slope 
was  constant  and  then  successively  relax  equality  until  the  maximum  was  reached. 
Clearly  the  unavailability  of  an  explicit  formula  carries  with  it  other  short¬ 


comings  . 
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